A FREE BOUNDARY PROBLEM ARISING FROM SEGREGATION OF 
POPULATIONS WITH HIGH COMPETITION 



Abstract. In this work, we show how to obtain a free boundary problem as the hmit of a fully 
non linear elliptic system of equations that models population segregation (Gause-Lotka-Volterra 
type). We study the regularity of the solutions. In particular, we prove Lipschitz regularity across 
the free boundary. The problem is motivated by the work done by Caffarelli, Karakhanyan and 
Fang-Hua Lin for the linear case. 



The problem we study in the present paper is motivated by the Gause-Lotka-Voltera model of 
extinction or coexistence of species that live in the same territory, can diffuse, and have high 
competition rates. 

Consider the equation, 



diffusion term '^-^ 
which models populations of different species in competition, where 

Ui{x,t) is the density of the population i at time t and position x; 

Ri is the intrinsic rate of growth of species i; 

di is the diffusion rate for species i; 

Oj is a positive number that characterizes the intraspecies competition for the species i; 

hij is a positive number that characterizes the interspecies competition between the species 
i and j. 

In the papers [22^ this model was studied initially without diffusion. These papers studied how 
species can survive or get extinct with time, depending on the interactions among them. Upon 
adding diffusion, Mimura, Ei and Fang proved that the existence of a stable solution depends 
on the shape of the domain and on the relations between the coefficients in the equation. The 
characterization for two species has been proved to be easier, while the three species interactions 
remain to be fully understood in these papers. 

In the sequence of papers by Dancer and Du [13l [HI [18] the authors decided to first understand 
better the steady case (time independent) in order to obtain results for the parabolic problem. In 
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1. Introduction 
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these papers, one can find sufficient and necessary conditions for the existence of positive solu- 
tions (ui, U2) and (ui, U2, u^) with explicit conditions on the coefficients Ri, Oj, bij for the following 
problem: 

2 



(1) 



with i = 1,2 and i 



-Aui = RiUi - aiUi^ 

Ui = on dQ, 
Ui > in 0, 
1,2,3, respectively. 



in f]. 



The spatial segregation obtained in the limit as bij — )• 00 of the competition-diffusion system was 
associated with a free boundary problem by Dancer, Hilhorst, Mimura, and Peletier in |1?J (i.e. 
in the case of high competition between the species). In llSj the existence and uniqueness of the 
solution to Problem ([T]) with just two populations has been studied using variational methods. 

Later in |7j, Conti, Terracini and Verzini proved that the limit problem is related with the optimal 
partition problem in A'^ dimensional domains. Since then, several papers by Conti, Felli, Terracini 
and Verzini [HI |9l [lOl [TTl [12] studied with a general formulation, the existence, uniqueness and 
regularity for the asymptotic limit of the following system, 

1 , 



-Au\ = f{u\] 



1"2 



in 



1"2 



in Vl, 



Ui 



on do., i = 1,2. 



where (t)i{x)4>j{x) = 0, for i ^ j. In these papers, the existence of a limit pair of functions (mi,M2) 
such that — )• (ui,U2) when e — )• is shown to have a tight connection with two different 
mathematical problems. Namely, 

a) to find the solution of a free boundary problem characterized by the conditions: 

-Aui = fiui)x{Ui>o} « = 1, 2, 

Ui{x) > inil, i = 1,2, 
ui{x)u2{x) = in $7, 
Ui = (pi on90, i = l,2, 

b) to find the solution for a optimal partition problem. 

The existence and uniqueness of solution for a type of free boundary problem of the form 

-Au = f{u)x{u>o} 
u{x) > in 17, 
u = on d^l 

with u bounded, was studied using variational methods by Dancer [16j. 
Then the regularity of solutions for the free boundary problem 



(2) 



Au' 



1 



-uj ^ Uj in 17, f = 1, . . . , d, 



u'>0 



in fi, i 



1, 



. d. 



^i(^) — 4'i{x) > on (9 17, i 
on517, i / j 



k(t>j 



,d, 
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was studied by Caffarelli, Karakhanyan and Lin in [3l H] with the viscosity approach. 

More specifically, in [1] the authors proved that the singular perturbed elliptic system Q has as 
limit, when e — t- 0, the following free boundary problem 

' Aui = when Ui > 0, i = 1, . . . , d, 

A(uj — Uj) < in rj, i = 1, . . . , d, 

Ui{x) > mfl,i = l,...,d, 
Ui Uj =0 in fi, i 7^ j, 
bi on (90, i = 1, 



(3) 



Ui 



,d. 



They also proved that the limit solutions Uj are Holder continuous and have linear growth from a 
free boundary point. Also that the set of interfaces {x : u(x) = 0} consists of two parts: a singular 
set of Hausdorff dimension n — 2; and a family of analytic surfaces, level surfaces of harmonic 
functions. 

The goal of this paper is to generalize the regularity results for the system ^ and ([s]) , presented 
in Sections 1 and 2 of [3], to the following nonlinear elliptic system of equations 



M^iul) = -ul^u'j, inn, i = l,...,d, 



ul > in ri, i = 1, . . . , d. 



hi (f)j 



on dn, i = 1, . . . ,d, 
on dil, i / j. 



where M denotes the extremal Pucci operator (see ([5j)), and to characterize the analogous limit 
problem Q for this case. We also address the existence of solutions of this system. 

We have chosen this problem, besides its intrinsic mathematical interest, in order to study a model 
that takes into account diffusion with preferential directions, so we are able to model situations 
with maximal diffusion. The choice of the operator is also related with its natural comparison with 
a non-divergence linear operator with measurable coefficients. 

The paper is organized as follows: in Section [2] we present the main theorems as well as some 
definitions and background. The proofs of these results are presented in different sections for the 
reader's convenience. Thus, Section [3] is dedicate to the existence of solutions u*^. The proof of 
Holder regularity up to the boundary for an equation of the type A^~(u) = f{x) with Holder 
boundary values in Lipschitz domain is a subsection of this section. Then, in Section |4j we prove 
Holder regularity uniform in e for u*^. The characterization of the limit problem as a free boundary 
problem has its proof in Section [5] The main result, the linear decay from the free boundary, is 
proved in Section |6j The appendices contain some technical proofs and results that could distract 
the reader from the essential ideas if they were to be presented in their respective sections. Thus, 
Appendix |A] has some elementary general properties of the Pucci operators. Appendix [B] has the 
Fabes and Strook inequality used in Section [4j Appendices [C} |D] and [E] contain the Alt- Caff arelli- 
Friedman Monotonicity formula and the proofs of some results used in the proof of Theorem 2.5 in 
Section [6l 



2. Main Results. Background and some definitions 



(4) 
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Let C be a bounded domain where d populations co-exist. Consider the following system of 
fully nonlinear elliptic equations with Dirichlet boundary data for 

M~{ul) = ^ul'^Uj, i = l,...,d, inO., 

^u\ = (t)i^ i = l,...,d, on dQ, 

where ul,{i = 1, . . . ,d) are non-negative functions defined in 0, that can be seen as a density of the 
population i, and the parameter ^ characterizes the level of competition between species. 

Each (pi is a non- negative Holder continuous function defined on dQ such that (f)i{x)(j)j{x) = for 
i 7^ j, meaning that they have disjoint supports. 

Here M~ denotes the extremal Pucci operator, defined as 

(5) M^{ijj):= inf aj^Aj (a;(x)) = A + A e,, 

ei<0 ei>0 

where .4.a,a is the set of symmetric n x n real matrices with eigenvalues in [A, A], for some fixed 
constants < A < A, and ej are the eigenvalues of the matrix D'^uj{x). 

We assume that u] are bounded, < n- < A^, for all i. Note that AAo; > 7W~(a;), thus tt| are 
subharmonic, for all i. 

Remark 2.1. Observe that if u is continuous and subharmonic in the viscosity sense then u is 
subharmonic in the distributional sense, meaning that 



/ An (/) da; := / u A(/> da; > V(/> > 0, G 
Jn Jn 



Our results in this paper are the following: 

Theorem 2.2 (Existence). Let e > constant, and Q be a Lipschitz domain. Let (pi be a 
non-negative Holder continuous functions defined on d^l. Then there exist continuous functions 
{ui, ■ ■ ■ ,u^) depending on the parameter e such that uf is a viscosity solution of 

M-~ (uj) = -u'j^'^^Uj, for all i = l,...,d. 



Theorem 2.3 (Regularity of solutions). Let e and (pi be as in Theorem 2.2 Let = (uf, • • • ,ii^) 
be solutions of Problem in i?i(0). Then there exist a constant a, < a < 1, such that for any 



£ (C"(Si(0)))'' and 



d < C{N), 



with N 



sup,- 



L°°(Bi(0)) 



and C{N) independent of e. 



In the limit as e — )• 0, this model forces the populations to segregate, meaning that in the limit the 
supports of the functions are disjoint and 

^ fi in the sense of measures, when e — t- 0. 

e 

The measure /i has support on the free boundary. Recall that the support of a measure /i is the 
complementary of the set {E : E the biggest open set such that ii{E) = 0}. 
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Theorem 2.4 (Characterization of the Hmit problem). Let (f)i be as in Theorem 2.2. If 
«G (C")°' is the limit of solutions of (Mj, then 



(1) M- (ui - Zk^. Uk) < 0; 

(2) (supp Ui)° n (supp {J2k^i ^fc)) = /or i = 1, . . . d; 

(3) A4^ (ui) = 0, when Ui{x) > 0, for x G 17 i = 1, . . . , d; 

(4) Ui{x) = (j)i{x), for X G dQ, i = 1, . . . , d. 

Theorem 2.5 (Lipschitz regularity for the free boundary problem). If u ^ (C"(i?i(0)))'^ is the 
limit of solutions of ^ in -Bi(O), and xq belongs to the set d{suppui) n Bi{0), then, without loss 
of generality, the growth of ui near the boundary of its support is controlled in a linear way and 
ui is Lipschitz. More precisely, there exist a universal constant C such that for any solution u, for 
any point xq on the free boundary: 

(1) supb^(^q)Ui <CR, 

(2) \\ui\\Lip{BR{xo)) ^ 

where C = C{n, ||w||^2(^j)) and R < \. 

The proofs are developed in the following sections. Although, the last three results are similar in 
spirit to the ones proved in [4J for the elliptic linear system of equations, our proofs use different 
techniques. 

2.1. Background and some definitions. In this part we present the definition of viscosity solu- 
tions, the spaces S(A, A, /) and S{X, A, /) and some results from the fully nonlinear elliptic theory 
(see f2] for the proofs and more detail). 

Here will denote the positive extremal Pucci operator, 

M^{uj) = sup aijDij{uj) = A + A Cj. 

^^-^^•A ei<0 ei>0 

Definition 2.6. Let f be a continuous function defined in Q and < A < A two constants. We 
denote by S{X, A, /) the space of continuous functions u defined in Q that are viscosity subsolutions 
of M.^{u) = f{x) in O, meaning that if xq ^ VL, A is a neighborhood of xq, and P a paraboloid 
(polynomial of degree 2) that touches u from above at xq, i.e. 

P{x) > u{x) X £ A and P{xq) = u{xo), 

then 

M+iP{xo)) > f{xo). 

In similar way, we denote by S{X, A, f) the space of continuous functions u defined in Vt that are 
viscosity supers olutions of Ai~{u) = f{x) in J7, meaning that if xq £ VL, A is a neighborhood of xq, 
and P a paraboloid (polynomial of degree 2) that touches u from below at xq, i.e. 

P{x) < u{x) X £ A and P{xo) = u{xo), 

then 

M-{Pixo))<fixo). 
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Remark 2.7. As in [2] we will denote by 5*(A, A, /) the set of viscosity solutions 

5(A,A,-|/|)n5(A,A,|/|). 

Definition 2.8. Let F : S x Q ^ M., where S is the space of real n x n symmetric matrices and 
n C M". We say that 

- F is a uniform elliptic operator if there are two positive constants A < A, called ellipticity 
constants, such that, for any M G S and x £ 

A ||iV|| < F{M + N)- F{M) < A ||A^|| 

for all N £ S non-negative definite matrix, where \\N\\ = sup|^|^| |A^3;| is the value of the 
maximum eigenvalue of N if N >0. 

- F is concave (convex) if it is concave (convex) as a function of M G S. 

Now, we recall the comparison principle for viscosity solutions, Corollary 3.7 in [2j, that states that 
a viscosity subsolution that is negative on the boundary has to remain negative in whole domain, 
and that a viscosity supersolution that is positive on the boundary has to remain positive in whole 
domain: 

Proposition 2.9. Assume that u G (7(17). Then, 

(1) u G ^(A, A, 0) and u <0 on imply u < in ^l. 

(2) u £ S{X, A, 0) and u > on imply u > in 

The following compactness result (Proposition 4.11 in [2j) follows from the closedness of the family 
of viscosity solutions of Problem ^ under the uniform convergence and the Ascoli-Arzela theo- 
rem. 

Proposition 2.10. Let {i*fc}fc>i be a sequence of uniformly elliptic operators with ellipticity con- 
stants A, A and let {uk}k>i C C(il) be viscosity solutions in Q of 

(6) Fk{D^Uk,x) = fix) 

Assume that F^ converges uniformly in compact sets of S x 0, to F, where S is the space of real 
symmetric matrices, and that is uniformly bounded in compact sets of Q.. Then there exist 
u G C(0) and a subsequence of {uk}k>i that converges uniformly to u in compact sets of $7. 
Moreover, F{D^u,x) = f{x) in the viscosity sense in Q. 

Below is the L^ Lemma that follows from Lemma 4.6 in [2j, using a standard covering argument. 
Note that it is enough to consider /+ instead of |/| due to the Alexandroff-Bakelman-Pucci estimate. 
Theorem 3.2 in [2]. 



Lemma 2.11. If u £ S{f^) in -Bi(O), u £ C(i?i(0)), f a is continuous and bounded function in 
Bi{0), and they satisfy: 

(1) hifij^(o) u{x) < 1 

(2) u{x) >0 in Bi{0) 

(3) ll/"^llL"(iJi(o)) ^ 

Then, if cq is sufficiently small, there exist d and e positive universal constants such that: 



{x £ Bi{0) : u{x) > t} 



< dt"^, for all t > 0. 
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Now, we recall the inequality that gives interior Holder regularity and that follows from the Harnack 
inequality for viscosity solutions: 

Proposition 2.12. Let uj G S{X,A, |/|) n ^(A, A, — |/|) with f a continuous and bounded function 
in Bi{0). Then, there exists a universal constant fi < 1 such that 

OSCb^(0)UJ < /iOSCijj(o)W + ||/||Ln(B^(o)) • 

The interior Holder regularity that we will use is a particular case of Theorem 7. 1 in [2] , and Sobolev 
embedding: 

Proposition 2.13. Let u be a bounded viscosity solution of Ai^{uj) = f{x) in i?i(0), with f a 
continuous bounded function in Bi (0) . Then there exists a positive constant C depending only on 
n, A, A such that to £ W'^'^{Bi (0)), for any p < oo, and so u G C^'°'{Bi (0)) for any a < 1, and we 

have 

ll^llci.<^(Bi (0)) - ^ ( II'^IIl°°(Bi(0)) + II/IIlp(Bi(0)) ) ■ 

Remark 2.14. 

(1) The same result under the same hypothesis is also valid for a general uniformly elliptic 
operator, concave or convex. 

(2) Observe that if f e then u G C^'". 



3. Existence (Theorem 2.2) 



To prove the existence theorem, Theorem 2.2, we will need a fixed point argument that can be 



found in [21], pg 280. We recall the result here for the sake of completeness: 

Proposition 3.1. Let a be a closed, convex subset of a Banach space B. Let T : a 
continuous function such that T{a) is a pre-compact set. Then T has a fixed point. 



a be a 



To apply the fixed point theorem, we need an existence result and regularity up to the boundary 
for a Bellman-type equation. In the following results, we denote by G the operator 



G[u;,] ■.= GiD^Ui,x) 



inf 

ast G C 



-UJi 



E 



M~ 



[UJi 



-UJi 



[ast] G ^A,A 

with Uj fixed positive continuous functions and ^a,a the set of symmetric nx n real matrices with 
eigenvalues in [A, A], for < A < A. The existence result is Theorem 17.18 in [21], but we also state 
it below in the adequate form for our purpose. 

Proposition 3.2. Let be a bounded domain in M" satisfying the exterior sphere condition for all 
X G dVL. Let Uj, j / i, be given functions, and Ogt symmetric matrices. Suppose that, for all s,t,j, 
there exists a positive constant fj, such that: 

Ost G C'^{Q) and \\ast\\c2(^n) < M-^/ 



G C72(0) and 



^u' 



C2(n) 



< ^A; 



and 
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0<A|Cp<a,t66<A|eP, and u]>0. 
Then, for any <pi G C{dQ), there exists a unique solution uji G C^(0) n C{Q,) of 

G{D'^u)i,x) = 0, in VL 
LOi = (j)i on dVl. 

We also need a generalization of the comparison principle that the reader can find on page 443 in 
|21j and comment on page 446: 



Proposition 3.3. Let n, w E C{^) fl C^(J7). If G[u\ > G[v] in Q and u < v in dO., then u < v in 



3.1. Holder regularity up to the boundary. The next Proposition is the Holder regularity up 
to the boundary for a viscosity solution of an equation of the type M.~{u) = f{x) in a Lipschitz 
domain. The proof we present here uses the comparison principle and an inductive construction 
of barriers. A different proof, by Luis Escauriaza, can be found in [TU], Lemma 3. Escauriaza uses 
the Fabes and Stroock inequality to estimate the Holder norm of the solution up to the boundary 
in terms of the L'^ norm of the right-hand side. 

This result is an improvement of Proposition 4.12 and 4.13 in [2] for general Lipschitz domains. 
Proposition 3.4. Let be a Lipschitz domain. Let uj G C'^{Q) n C{^) he a viscosity solution of 

{M~{ijj) = fix) in i}, 
Ld = (j) on dQ, 



(7) 



with f >0 and f £ C{Q) n C^{n), cj) G C^^idn). Then uj G C^iQ), where 7 = min(§,a) and a is 
the universal Holder exponent for the interior regularity. 

The proof of this Proposition follows the same lines of the proof of Proposition 4.13, in [2]. Once 
the interior regularity and the regularity for an arbitrary point in the boundary is guaranteed the 
proof is basically the interplay of these two results depending on how close two points are compared 
to the maximum of their distance to the boundary. The interior regularity comes from Proposition 
2.13 To prove the regularity for an arbitrary point on the boundary. Proposition 3.10[ we need the 



couple of Lemmas that follow. 

The first Lemma establishes the decay of the solution in concentric balls centered at an external 
point in the outside cone. The proof uses a standard comparison argument and the use of a barrier 
function. 

Lemma 3.5. Let 0, be an Lipschitz domain and C an external cone centered at xq G dQ, with 
some universal opening. Let y £ C be the center of the balls C B'^ C B^ such that B^ G C , and 
that dist(9i?^ H Q.,xq) > 5 > 0. Let u be a solution of M.~^{u) > in the viscosity sense in $7, such 
that u < 1 on B^ n Q and u < on B^ Ci dO,. Then, there exist A > such that 

u{x) < A < 1 in B'^nQ. 



Proof. Since the domain is Lipschitz there exists a cone C with opening equal to p, such that for 
any point of the boundary we can place the cone with opening p and vertex at that point such that 
cno = 0. 
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Figure 1. Reference decay. 



Without loss of generality, take the cone with origin at xq G and with axis Cn , 



c 



n-1 



and consider the origin as the center of the balls B^, B^, B^ as illustrated in Figure [T] 



Applying Lemma E.2 with M = ^ where r is the radius of B^ and ^ equal to the radius of B^ we 
obtain a supersolution -;/;, A^''"(^) < in B'^\B^, such that iplx) = in dB^ and ^{x) = 1 in dB^. 
And so 

u{x) < ipix) on d{B^ n n) 
Applying the comparison principle stated in Proposition |2.9| we can conclude that 

u{x) < ip{x) in B^ nfl 




Let A = ip{x) < 1 for x G dB^ CiQ. As ^p is an increasing function, see Figure [2| we can conclude 
that for x & B'^ CiQ, we have that u{x) < A < 1 as we claim. □ 
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The next Lemma is important for the iteration construction. Basically we prove that if the boundary 
data is bounded from above in half of the unit ball centered on a boundary point, and the function 
is bounded in the unit ball then in one fourth of the ball, the function decays by a fixed value. 

Lemma 3.6. Let Q be a Lipschitz domain and C an external cone with universal opening. Let 
S dVt he the origin of the cone. Let v he a solution of A4~^{u) > in the viscosity sense in fl, 
such that V <l on Bi{0) D n, v{0) = and v{x) < (l)^ on Bi{0) D dn, for some /3 > 0. Then, 

there exists a constant (^)'^<^<1 such that 

v{x) < n in Bi{o)nn. 



Proof. By hypothesis, 



Let 



v{x) < ( 2 ) ' ^ ^ dnnBi{o). 



uj{x) 



y(=^) _ (1)13 

for xGSi(0)nJ7i, 



where fii = {2x : x G ft}, uj satisfies all the hypotheses of Lemma 
(1) xo = 0; 



3.5 



with: 



(2) C the uniform external cone with axis without loss of generality equals to e„ axis; 

(3) B^ = Br{y) withy = (O,--- ,0,-^^5) and r < dist(2/, aC), r = ^; 

(4) i?2 = i?_^(y); 

(5) B^ = B,{y). 

Observe that B'-^ C Si(0) and that -Bi(O) C B"^ . Then by Lemma 

uj{x) < A < 1 in n r2i 



3.5 



and so we have also that 



^^^^^ <\^v(^)<\{l 



(5)^^ ^ M — 1- Then, we obtain that 



/9 



for xGSi(0)nni, 

2 2 



v{z)<ii z£Bi{o)nn. 



Now we will be able to prove an iterative decay illustrated in Figure [3) 



□ 



Lemma 3.7. Let il. he a Lipschitz domain and C an external cone with universal opening. Let 
G do, he the origin of the cone. Let v he a solution of M.~^{u) > in the viscosity sense in O, 
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Proof. By scaling and dilation, define for x £ Bi (0) 

f^k - (4F2) 



Since oj satisfies the hypotheses of Lemma 3.6, considering this dilation and scaling, we see that 

^(x) < A < 1 in Bi{0)nn 
and so, like in the previous proof, we have also that 

' _ T_H <A^^(i^)<V. + (i-A)(^^J for xGB.{0)nn. 



Therefore, 



v{y) < Hk+i for y e B^ (0) n n, 

4^+1 



with {^^2)^ — f^k+i < fJ-k- This finishes the proof. □ 

Remark 3.8. Observe that the decay at each step of this iteration is constant and equal to Co (|^)" 
for a much smaller than f3 and Cq a large positive constant. In fact, there exit constants Cq and 
a such that jJ-k ^ Cq {^^^ . By induction, for k = the result is true by Lemma 3.6 for Cq > 1. 
Assuming the result valid for a general k, we have that 

^,k+l = A^. + (1 - A) f -l-V < — i- f 4" ACo + (1 - A) ^ 
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Take a = e/3 and such that 4"A < 1 then 



for Cq large constant since, 



la 1 

4 ACn + (1 — A) — 7^ — TTT^ r = 4 A Co + (1 — A) — Trr, — — rTTT, T ^ Cq 

^ ' 2/^-2o4A:(/3-a) ^ ^ ' 2/3(l-2e)4A:/g(l-e) — 

<1 



(1 - 4"A)2/^(i-2^)4'=/5(i-^) ■ 

Remark 3.9. Note that Lemmas 3.5 - 3.1 are valid for v a viscosity solution of problem 

Proposition 3.10 (Holder regularity up to the boundary). Let he an Lipschitz domain and C 
an external cone with universal opening that only depends on the domain. Let xq G be the 
origin of the cone. Let v be a viscosity solution of problem ^ such that \v{x)\ < 1 on Bi{xo) n Q, 
v{xq) = and \v{x)\ < \x\^ on Bi{xq) n dCl. Then, there exist constants C > and a << /3 such 
that 

andC = Ci\\f\\L^). 

Remark 3.11. The constant C in the previous proposition would depend on e if this result was 
to be applied to our main problem. The uniform Holder regularity in e will be proved in the next 
section and does not depend on this result though. 



Proof. Assume by translation invariance that xq = 0. Note that v{0) = v{xq) = 0. So if we prove 
that 

\v{x)\ < C \x\°' , 

the result follows. Observe that on the boundary the regularity of the boundary data gives the 
result directly. Let k be such that 

1 , , 1 

< T < 

4fe+i - I I - 4k- 



Using Lemma 3.6 followed by Lemma |3.7| (see Remark 3.9) we can assume that for x such that 
|x| = yO < ^ we have: 

v{x) < flk- 

Then, taking in account the previous remark, we also have, 

1 



v{x) < Co 



4k 



but then 



vix) < Co4° 



4^+1 



< Co4"|xr. 



To obtain the other inequality observe that —v{x) < 1. So if we consider instead of —v the function 



u{x) 



2n\ 



\x\ — v{x). 
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we have that, for x E -Bi(O), 
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2nX 



Observe that we have as weU, for x G 317 n -Bi(O) that 



ijj{x) 



2nX 



(x) < 



2n\ 



+ C\\xf. 



Therefore, we can apply the comparison principle for uj and a barrier function as in Lemma 3.6 



Repeating the same construction as in the proof of Lemma 3.6 and Lemma 3.7 we obtain also that 
and so, in an analogous way, we obtain that 

-v{x)<c\pr. 

with C depending of ||/||^oo. This completes the proof. □ 



With this result that is the analogous of Proposition 4.12 in [2], the proof of Proposition 3.4 follows 
as the proof of Propositon 4.13 in [2]. 



3.2. Proof of Theorem 2.2 , In this section we finally present the proof of Theorem 2.2 stated in 
Section [2j 

Proof. Let B be the Banach space of bounded continuous vector-valued functions defined on a 
domain Vt with the norm 



= "l'^^ ( sup \Ui{x)\ ] 



Let a be the subset of bounded continuous functions that satisfy prescribed boundary data, and 
are bounded from above and from below as is stated below: 

a = {{u\,u\, ■ ■ ■ , li^) : uf is continuous, = (j)i{x) when x £ dQ, < ul{x) < sup ||(/)j||^oo} 

i 

0" is a closed and convex subset of B. Let T be the operator that is defined in the following way: 
7" ((-"!, -"Ir • • ,Ud)) = {vhv2,-- ■ if K,""!,- • • ,u^d) {vi,V2,-- ■ ,Vd) are such that, 

1 
e 



M-{vl) = -Yl 
i = 1 



l,...,d 



in J7 



, . . . , ^, 



, d, in dQ, 

in the viscosity sense, where j / i are fixed. Let g(u'^) = ^ X^jyi ^-^d observe that each of 
the previous equations of the vector u*^ = • • • , u^) take the form: 

M-{vt) = vlg{vi^). 

Observe that if T has a fixed point, then 

T u|, • • • , u%)) = {u\,ul, ■■■ ,ul) 

meaning that = (pi on the boundary and that 



in Q, which proves the existence as desired. 



14 VERONICA QUITALO 

So in order for T to have a fixed point we need to prove that it satisfies the hypotheses of Proposition 

(1) T{a) C a : We need to prove that there exists a regular solution for each one of the equations 

iM-{vl)-vlg{M') = Q, k = l,...,d, in n 
\vl = (f>k, k = l,...,d, in dft. 

Observe that if such a regular solution exists, the comparison principle is valid by Proposi- 
tion [ 



To use Proposition 3.2 we can rewrite the differential equation in the form 



F{vl)= inf (a,,AA-^fc9(u^)) = 0, k = 1, 

aij G Q 



since, by density, taking the infimum in Q is equal to taking the infimum over M. Let 
be a family of smooth domains contained in such that as / — )• oo. Since we do 

not have the desired regularity on the coefficients of v^. we will first consider the boundary 
value problem in each smooth domain $7; with the following regularized equation, 

F\vI)= inf {a^jDijvl-vl{g{u')*ps)) = 0, k = l,...,d, 



aij G 
[aij] G A 



A,A 



where ps is an approximation of the identity. And we will consider suitable boundary data 
that converges to the original boundary data when we approach the original domain. 



So now, by Proposition 3.2, there exist a unique solution {vl)f £ C^{n) n C°(J]) with 
K)f = (<^fc)« on dCli, k = l,---,d. 

Taking the limit in / we obtain the existence of {vlY, solutions in ft with boundary data 
equals to (p^- From Proposition 3.4 we also have that there exists a universal exponent 7 
such that (v^.)^ is Holder continuous up to the boundary, (f^)^ G C'^(rj). 

As we have ^ F uniformly due to the properties of identity approximation, we can 
conclude that 

F'ivf) ^ F{vl) 



by Proposition 2.10, this is, that G ^^{fl) is the solution for our problem F(vl.) = in 
fl, for all k. 

We need also to prove that for each A:, < < supj ||(/)j||^oo . 

Suppose, by contradiction that there exists xq such that vl{xo) < 0. Since v^. = (pk on dfl 
and 4>k is non-negative, xq must be an interior point. But, if has an interior minimum 
attained at xq then there exists a paraboloid P touching v^. from above at xq and such that 

M-{P{xo))>0. 

Since, (?(u'^) > 0, on the other hand we have that 

M-{vUxo)) = vl{xo)gin'{xo))<0, 

and so we have a contradiction. 
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Analogously, suppose, by contradiction that there exists xq such that f^(xo) > supj ||(/).j||^oo . 
Then, by the same reason as before xq must be an interior point and 

M~{vl{xo))<0. 

Since, g{u^) > 0, we have that 

M-{vl{xo)) = vl{xo)g{u'{xo))>0, 
and so we have a contradiction. 

T is continuous: Let us assume that for each fixed e we have ((nf )„, • • • , — )■ (uf , • • • , 
in [C(r2)]'^ meaning that when n tends to oo, 

max ||(n,^)„ -M^ll^^ 0. 

l<i<a 

We need to prove that 

||r((nf)n,--- ,(M^)n) - T(nf,-- - ,u^)||[c(n)]<i ^ 
when n — )• oo. Since, 

T (K)n,, • • • , K)n) = {{vDn, ■■■ , (w^)n) 

if we prove that there exists a constant C, independent of i, so that we have the estimate 



< C max 
j 



L°° 



m)n-vi\\^. 

the result follows. For all x G $7, let a;„ be the function 

W„(x) = {V^i)n{x) -Vlix), 

and suppose by contradiction that there exists y £ ^ such that 
Wn(y) > r^Kmax (np„ - Uj 

for some large K > 0, where r is such that ^1 C Br{0)- We want to prove that this is 
impossible if K is sufficiently large. Let hn be the concave radially symmetric function, 

hn{x) = 7(r2 - 

with 7 = K mayii — ii^||^oo- Observe that: 

(a) hn{x) = on 55,(0); 



(b) < K maxj (np 



for all X in Br{0); 



(c) = Unix) < hn{x) ou 517, since (w|)n and v| are solutions with the same boundary 
data. 

So, since we are assuming ([s]), there exists a negative minimum of hn — Wn- Let xq be the 
point where the negative minimum value of hn —Un is attained, hn{xQ) — a;„(xo) < 0. Then 
we have that for any matrix A £ ^a,A) 

OijDij {{hn - ujn) (xo)) >0 and aijDijhn{x) = -270^ < 0. 
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Moreover, 

M + {LOn) > M-{{vl)n{x))-M-{vl) 



iiv!)n-vl) 



ek)"-<eH-(-P")) 

> \{{{vt)n-vl)Y,{'^%-vl{d-l 



adding and substrating -w,- ^j-ti{u'j)n- Then, if afj are the coefficients associated to uj^, i.e. 



then 

< a^" Ai (^n - w„) (xo) 



< -2 7C(^o) - I ( m)n - vf) (xo) 5]^-)„(xo) - vl{xo){d - 1) ||(np„ - u^^\\^^ j 
because < hn{xo) < a;„(xo) and Ylj:^ii'^j)n{xo) > and so 

-\{iv'r)n-vl){xo)Y,iu^^Uxo)<0. 

Taking K > supj H^iH^oo , we obtain that 



< a^/Aj(/ln-C^n)(xo) <0 



which is a contradition. 



(3) T{a) is precompact: This a consequence of (1), since the solutions to the equation are 
Holder continuous, and this set is precompact in a. 



This concludes the proof. 



□ 



4. Regularity of solutions (Theorem 2.3) 



Following the proof presented in [1] , we will prove the decay of the oscillation of the solution under 
certain particular hypotheses and then use this result to prove the uniform C° regularity for a 
solution of the system of equations. 

4.1. Some Lemmas. We first establish some conditions under which the maximum of a function 
in a smaller ball is lower. 

Lemma 4.1. Letu^ he a solution of Problem . Let Mi = max^-g^^^o) (•^) o^^dOi = osCj,.gBj(o)'^i(^) 
If for some positive constant 70 one of the f (mowing hypotheses is verified 



(1) {xGSi(0):n^(x)<M, -7oOi} 

4 

(2) {xG Bi(0) :7W-K^(x)) > 70 Oi} 



> 70 
> 70 
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(3) {x£Bi{0):M^iu^ix))>-foul{x)} 



> 70 



then there exist a small positive constant cq = co(7o) such that the following decay estimate is valid: 



max uUx) < Mi — Co Oi 
xe-Bi(o) 

3 



and so 



oscb^(0)uI{x) < Co Oj, 



with Co < 1. 



Proof. 

(1) By contradiction, assume that for all co small, exist a point xq € Bi{0) such that 

4 

ulixo) > Mi - Co Oi 

and let 



Vi{x) 



Mi-ul{x) 



Oi 



Vi satisfy the following properties, with fi[x) 



(a) inf^^ 



(0) 



CO 



<mfB.(0)^<l 



(b) >[-(|)</.(x) 

(c) Vi{x) > in Bi{0) 



Then we can apply the L'^— Lemma stated as Lemma 2.11, to conclude that 



{xeBi(O) : ^ >t} 



Co 



< dt~ 



for d, 6 universal constants and for all t > 0. If t = — then we have 

' Co 

{xeBi{0):Vi{x)>jo} <d{^)-^ 

4 Co 

But, taking in account the hypothesis we have: 

70 < 



{xGB^{0)■.v,{x)>-fo} <d{^)~^ 

4 Co 



S+1 

Since cq is arbitrary, if Cg < we have the contradiction. 



(2) If u- is a solution for M (n-) = f{ul) there exists a symmetric matrix with coefficients 
aij{x) with 

x\^\^ < aijCiCj < m"" 

such that 

aki{x)Dkiul(x) = f{ul{x)) 
For that particular matrix consider the linear problem with measurable coefficients in i?i(0) : 

La {v{x)) = g{x) 

and let G{x, •) be the respective Green function on i?i(0) such that for all x G i3i(0), 



v{x) 



Bi(0) 



G{x, y)g{y)dy + boundary terms 
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(see page 415 in [H]). We conclude 



Mi-ul{x) > 



BiiO) 



G(x,y)(-/K))dy> / (G(x,y))(/K))dy, 



A, 



for Ai := {x € -Bi (0) : M~{ul{x)) > 70 Oi}, since 

4 

aki{x)Dki{Mi - u\[x)) = -f{u\{x)), 
and the boundary values are positive and G{z^ y)fiul{y)) > for all y. Since by hypothesis 



{xGBi{0):M-{ul{x))>-foO,} 



> 70 



and due to Fabes-Strook inequality, (see Lemma p.l and Theorem 2 in 120] for more details) 

I A; I 



G{z, yho Ojdy > 70 Oi 



A, 



We claim that: 



Si(0) 



4 



G{z,y)dy. 



\f G{z,y)2nMy>C. 
is 1(0) 



2nA 

So again by hypothesis and due to the claim (that we will prove later) we have 



70 Oi 



\Ai 



Bi(0) 

3 



SO 



i?i(0) 

4 

Mi-U\{x) > -foO^ 

with Co < 1. 

To prove claim ^ we argue that 



G{z, y)dj > 70 Oi 



70 



5i(0) 



c 



70 



i?i(0) 



B 



Giz,y)2nAdy> [ G{z, y)2 ly^an] dy = ^ 
1(0) Jbao) \i J 42 



Ul > c 



for 2; interior, since aijDij{^ 



-2 Y2i (^ii s-iid — 2nA < —2^ - Oji < —2n\ 



(3) Let 



and 



Ai = \xe 5i(0) : 7W-(uKx)) > 70 -"^(3^)} 



Hi = {xeAi: ul{x) < 



Mi 



and consider the two possible cases: (a) \Ai\Hi\ > | \Ai\ and (6) |A\-f^i| < 5 |A| 
(a) If |A\-f^i| > 5 |A| then as 

X G A\^^ : A^-«(x)) > 70 ^1 C |x G 5i(0) : M-{u^{x)) > 
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since < we can conclude that 



X e Bi(0) : M~{u\{x)) > 70 



Then we have the decay by ^ with 70 replaced by 
(b) If \Ai\Hi\ < l\Ai\ then as 

Hi = !^xeAi: ulix) < ^1 C {x G Bi(0) : ut{x) < - /3oO,} 



for f^o < i§- and 



m = \A,\ {A,\H,)\ = \A\ - \{A,\H,)\ > ^ > y 



we have 

{x G 5i (0) : ulix) <Mi- ToOi} 
for To = min(/3o, 7o)- The decay follows by (1). 



> 



To 



□ 



Next Lemma states that if all the oscillations are tiny compared to just one that remains big, then 
the largest oscillation has to decay due to an increase of the minimum in a smaller ball. 

Lemma 4.2. Let v!^ he a solution of Problem ^ in i?i(0). Let 

0} = OSC^gB^(o)M^(x). 

Assume that for some 5 > 0, sufficiently small, 

then 0\ must decay in Bi{0), that is, there exist fi < 1 such that 



Of < nO\. 




Figure 4. Illustration of our hypotheses about oscillations. 

Proof. Let u be the solution of the problem 

M-{oj{x)) = 0, x£Bi{0) 
u;{x)=ul{x), xGdBi{0) 
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Since uf — a; is a subsolution for the positive Pucci extremal operator 

M+{ul-uj) > M~{ul) + M+{-uj) > 0, 
and nf — cj + ~ ''^i) ^ supersolution for the negative Pucci extremal operator, 



< M-{u\) + Y,M+{Mi-v^^ 

< M-{u\) -^M-{u\) 

by the maximum principle for viscosity solutions and our hypotheses we get, 
u\{x) <uj{x)< u\{x) + ^(Mj - ulix)) < u\ + 60\. 

Since w G 5*(A, A, 0), uj decays, namely 

osc^eBi(o)^(a^) < osc^gs^(o)a;(3;) 

2 

and from this inequalities, since 

max uj — min cj < max u\ + 50\ — min oj < max u\ + 60l — min u\ , 
we can conclude that 

oscB^(o)a; < oscb,(o)^ + ^Ol- 

In a analogous way. 
So, then, 

oscbi(o)^ < oscBj(o)a; + (50{ < fi osc^^Bi{o)'^{x) + SOl < fi {oscb^(o)u\ + 50l) + 50l. 

1 2 

Simplifying, 

oscB^(o)nf < {ii{l + 5) + 6)0\, 

2 

which concludes the proof by taking 8 sufficiently small. □ 



4.2. Proof of Theorem 2.3, In this section we finally present the proof of Theorem 2.3 stated in 
section [2l 

Proof. We prove this theorem iteratively. We will prove that the oscillation of u*^ will decay, by 
some constant factor /i < 1, independent of e when it goes from -Bi(O) to B\{Q) for some A < 1 
also independent of e. Meaning that, we will prove that there exist two constants < A,/i < 1, 
independent of e, such that for alH = 1, . . . , d, we have 

OSC^&Bx(0)U'i{x) < floSC^^Bi(0)Ul{x). 

The regularity of each function will follow from here in a standard way using Lemma 8.23, in 
|21j . Since, what matters is the ratio between oscillations 

OSCj;^(0)^^(x) ^ ^ 
OSCb,(o)<(x) ~ 
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the result will hold true if we prove this decay for the normalized functions, which satisfy the same 
equation with a different value of e. Thus, consider lif , u|, • • • , solutions of Problem Q on i?i(0) 
and the renormalized functions u\, - ■ ■ , n^, 

ul{x) = pul{x), X e Bi{0), i = l, ■■■,(!, 

with p = er^- These functions are bounded from above by one and satisfy 



M-{ul{x))= ul{x)Y,ut{x), i = l,---,d. 



4^ i^l 
e 

Briefly, the iterative process consists of the following. We prove that in Bi (0) at least the largest 

4 

oscillation decays. Without loss of generality consider u\ the function with the largest oscillation. 
Then, there exists p < 1 such that 

OSC^(zB^(0)Ul{x) < 7xOSC^eB,(o)Uf(x). 

4 

Then, we consider the renormalization by the dilation in x : 

ul{x) = p ul (^xj , xeBi{0), i = l, ■■■,(!, 



with 



1 

P= ^>1- 

max^gs^(o),fc=i,.,.,d ul{x) 



Observe that these functions are solutions of the system 



So, basically, they are the solutions of an equivalent system with a different e, still defined on -Bi(O). 

We start all over to prove that we have again the reduction of the next largest oscillation, when we 
are in Si (0). We call the new function with the largest oscillation u\. So we have that there exists 

4 

p < 1, independent of e, such that 

osCa.gBj(o)<(x) < posc^gB^(o)M-(3^) =^ osc^gB ^ {o)ul{x) <jiosc^(,B^{o)ut{x). 

4 J2 4 

Considering what has happened in the previous step we can have several scenarios. Either we 
have the reduction of the oscillation from -Bi(O) to B i (0) of just one of the functions, or two or 

4? 

more. (Later, when there is no possible confusion, the function with the largest oscillation at each 
iteration will always be denoted by uf). 

We repeat this process taking the renormalizations 

ul{x) = ^ u\ (]^x\ , xESi(0), i = l,---,d, 

until eventually we have the reduction of all the oscillations, or until the largest oscillation is 
much larger then the other oscillations. If that is the case, then eventually we will have that, for 
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Figure 5. Decay iteration. After the renormalization the oscillation of the first 
function decays while the others remain the same. In the original configuration we 
register that decay and we proceed with the next renormalization. 



Oj = OSCBj^u'jix), 



In that case, using the Lemma 4.2 , Oi must also decay and we have the result. Thus, after repeating 
this iterative process a finite number of times we obtain that for some A < 1, /i < 1 and every 
i = 1, . . . 

For simplicity, we will still refer to the renormalized functions bounded by one, as n^, and also _ 
will still be denoted by e in each new step. Although, e can be bigger than the one in the first 
steps, depending on the number of steps needed, it will eventually remain smaller than one since 
the renormalization after each dilation will be a multiplication by a factor smaller than one. 

In what follows n^-^ denotes at each renormalization, the function that achieves the maximum value 
1 and u\ is the function that has maximum oscillation. Naturally, they can be or not the same 
function. 

Below follows the proof of decay for the renormalized functions in all possible cases. 
Case 1: Let us assume that ^ > 1- 

We are going to use the following argument: observe that, if there exists k such that: 



{x G 5i : ulix) > 70} 



> 70, 
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then since for any j ^ k 
we have, 



M-{u]{x)) = -yAx){ul + ---) 



> 70- 



And so by (3) in Lemma 4.1 we can conclude that 0^, for all j k decays. Let < 7 < ^ be a 
fixed constant. 



(1) If max^,. uf > 7 and Oi does not decay then by (1) in Lemma 4.1 we can conclude that there 
exists 70 a positive small constant such that 



|x G Bi : ti|(x) > maxn| — 70O1 ^ ^| 



> 1 - 70 > 



2' 



so then as we observed before, we have for all j ^ 1 that, 



4 -' 2 J 



> 



7 



And so by (3) in Lemma 4.1 we can conclude that Oj, for all j ^ 1 decays. 



(2) If max^ < 7 then, since Oi < max^^ < 7, all oscillations are smaller than 7 (see Figure 





Figure 6. All oscillations are smaller than 7. 

Then in particular the function Uj^ has oscillation smaller than 7. Thus, 

u%ix) > 1 - 7 
and so for all j 7^ jo, we have for all j 7^ jo, 



M- 



(np>(l-7H. 



4.1 



we can conclude that for all j / jo, Oj decays. In particular. 



So again by (3) in Lemma 
the largest oscillation decayed 

Case 2: Let us assume that ^ < 1- Let = ];■ 

Observe that, since all the functions are bounded from above by one and are positive, we have that 
for all i ^ io, and io an arbitrarily fixed indix, 

(10) eu%ui<M-{uf)<eui{d-i). 

In a more general point of view we have for any i = 1, . . . ,d, 

(11) 0<M~{ul)<d, 

We will use one or another expression according to convenience. There are two possible cases: 
either (1) | < Oi < 1 or (2) d < ^. We prove that in the first case all the big oscillations decay. 
So after a finite number of steps, (and if all the functions did not decay yet in the mean time) 
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we will have that all oscillations are less than |. In this case, since the function that attains the 
maximum will have also oscillation less than |, we can prove that all the other functions decay. In 
this way, we will eventually be again in the case of 

and as before we have the result. The proof of the decay in each case follows. 

(1) Assume that t < Oi < 1. There exists an interior point y G Bi{0), such that 

4 

m.inu\ H — < u\{y) < maxuf . 



By (11) we conclude that the equation for ul has right hand side continuous and bounded. 
Then, by regularity, Proposition |2 . 1 3 we can conclude that there exists a universal constant 
C such that for all x G (y) 

16 C 

K{x)-ul{y)\ < llV-ufll^^ ^ < ||nf||pi,„ ^ < ^. 
And so, for all x G B^_(y) we have that 

16 c 

u\{x) < maxn| . 

16 



Then by (1) of Lemma 4.1 , we can conclude that 

osc^^(o)^ ^ /^osc5j(o)^ with fi < I. 

With this argument we can conclude that all the functions with oscillations bigger than | 
decay. Moreover, repeating the same argument we can conclude that all the functions with 
oscillations bigger than | decay. 

(2) Assume that Oi < Observe that in this case Uj^{x) > | for all x € Bi{0) 
Then for all i ^ jo the function satisfies the equation 

^eul<M-{ul)<9ul{d-l) 

Meaning that, 

M-iui) ~ eul 

Consider Vi, i ^ jo, the renormalized functions u^, by the renormalization 

Vi{x) = ^ — ul {x) . 

max^gBj(o) ul{x) 

Observe that Vi has maximum 1 and satisfies in Bi{0), 

9^v^ < M~ivi) < evi{d-l). 

For each function Vi, i ^ jo, we can have two situations. Either (a) Vi{x) > \ for all x or (b) 
Vi{x) < ^ for some x in Bi{0). In both cases we will prove that Vi, i ^ jo decay in Bi{0). 



(a) If Vi{x) > i for all x G Bi{0). Then 



Oi = osc^.eij^(o)Wi(x) < ^. 
We claim that: there exists a universal constant N such that 6 < NOi. 
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Assuming that the claim is true, observe that 

< M~{vi) < dNOi, 

which imphes that we can consider the function uj with oscillation 1 in -Bi(O) defined 
by 

Vi{x) - minxes (0) Vi{z) 

= • 

and that satisfies in -Bi(O) 

< A^- (w) < dN. 



By regularity, Proposition 2.13 there exists a universal constant C depending on A'^ 
and d such that, 

|a;(x) - uj{y)\ < \\uj\\fji,c \x - y\ < C{N, d) \x - y\ . 
Assuming that u didn't decay, let y G ^i(O) be such that uj{y) < ^. Then, if 
a = min (^dist(y, dBi), jcW:^)) ' ^ ^ ^^(v) 

\uj{x) - u}{y)\ < J 

And so, for x G B^{y), 

,,113 
u:{x)<- + -<-. 

Then by (1) of Lemma 4.1 , we can conclude that for all i ^ jo, there exist fi < 1 

OSCB^{0)Vi{x) < fiOSCB^(0)Vi{^) oscb^(o)uI{x) < ^ osc^^ (o)n^ (x) (x) . 

Proof of the claim: consider by contradiction that 6 > 2n\ 24 Oj. So, for all x £ -Bi(O), 

3 3 

-Vi{x)e > -9> 2nX9 0i. 

4 8 

Observe that, if m = iiimz^Bi{o) "Viiz), m < t < m + oscBi(o)Vi{x) is a positive constant 
to be chosen later and 

P[x) = 8 0^ \xf + t 

we have that for any x, and any t, 

M-{vi) > ^Vi{x)e > 2nX9 0^ > M-{P{x)) = 2nX80i. 

and so P can not touch Vi from above at any y, since it would contradict that the Vi 
is a subsolution. But since for t = m and |x| = 1, 

P{x) >m + SOi and P(0) = m, 

P crosses the function Vi, and so it is possible to find t such that P would touch the 
function Vi from above at some point, which is a contradiction. 

(b) If Vi{y) < 2 for some interior point y of Bi (0) we proceed as before and use regularity. 

4 

If the function never attains a value less than ^ in Bi (0) then the function has decayed. 

□ 



With this uniform bound in the Banach space one can conclude that these sequence converges 
uniformly (up to a subsequence) to a vector of functions u. 



26 veronica quitalo 

5. Characterization of limit problem: Proof of Theorem 12.41 



In this section we will_assume without loss of generahty that A = 1 . Observe that if u*^ is a viscosity 
solution o: 
such that 



solution of Problem (|4|) then there exists a subsequence still indexed by e and function u G (C")*^ 



u*" — 7- u uniformly. 

The following Lemma characterizes the Laplacian of the limit solution. 

Lemma 5.1. If uG (C°)°' is the limit of a solution of ^ then Aui are positive measures. 

Proof. Let <j) be positive a test function. Then 

0< fM-{ul)4>< fAul4>= [ A(i)ul^ ! A(t)Ui= [ (pAui 



and so Am is a positive distribution which implies that it is a nonnegative Radon measure. □ 



Now we are ready to prove Theorem 2.4 (stated in Section [2j). 
Proof. (1) Observe that 

^ M~{ul)+M+l-Y,ul] <M-{nl)-M- 
< M-{ul)-Y,M-{ul)<^. 

k^i 

As u\ — Ylikjf^i'^k ~^ — X^fc-^jUfc when e — )• uniformly and 5(A, A,0) is closed under uniform 
convergence, Ui — Ylk^ti'^k is a supersolution of A^^: 

\ k^i I 

(2) If Ui(xo) = ao > for any i, i = 1, . . . , d then for 5 < ^ there exists an eo such that for e < eo, 



2 

a) \ul{y) - uI{xq)\ < ^ in B2h{xo)] 



^ < — 6 < uI{xq) < + 6 < So by Holder continuity there exist h > such that: 



b) uj{y) > ^ in a ball of radius 2h and center xq. 
Observe that applying Green's Identity with a function u and the fundamental solution 



na;„(2-n) Vlx-xol"-' 12/11"" V 
we obtain the inequality 

h^-f Audx<C'f {u{x) -u{xo))dS, 
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where C is a constant just depending on n. Prom the equation for uj, we obtain 

< / A<dx<-g/ K(y)-<(xo))dS 



< 



'Bhixo) "' JdB2h{xo) 

C ao 
47^' 



Thus, 



which imphes that, 



/ 

•I Bhixo 



) 2 



4^0 



when e — 7- . By subharmonicity, 

and so ^^.^j ui.{xq) = and this proves the result. 

(3) To prove that Ai~{ui) = 0, when Ui > assume the set up in the beginning of (2) for Ui{xo). 

We need to prove that — — ^ uniformly when e — >■ in order to use the closedness of 
^(A, A). Since, 

A (^5^^ > M- > J ^ 0' 

^ is subharmonic. So, if we prove that ^ ^ in L^(5/i(a;o)), then for y G -B/i_^/(xo), 

5Mi^</ ^^dx^O, 



e 



and so ^^^f^ ^ convergences uniformly in a compact set contained in Bh{xo). Recall that we 
proved in (2) that J^fe^jU^ — >■ uniformly in a compact set contained in Bh{xo) and we have that 



E ^ ^ ^ E ^ ill the sense of distributions. 

So, since 

— ^ < =M {ul)<2_^M u% 



we conclude that 

^^^0 inL\ 



^k I ' 
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Then, as we said, we have that 

> uniformly in a compact set contained in Bhixo). 

As uf — > Ui uniformly and are bounded, we finally conclude that 

u| > uniformly in a compact set contained in Bh{xQ). 

Proceeding analogously with Uk, k = 1, • • • ,n we conclude that the limit problem is 

M~{ui) = 0, Ui > i = l,...,d. 

(4) To prove the last statement we will construct an upper and lower barrier. 

(a) Consider as upper barriers the solutions of the d problems, {i = 1, . . . ,d): 

M~{u*) = mCl and u*{x) = (l)i{x)x{,pi{x)^o} ^^90. 
So we have that for alH = 1, . . . , d, 

M-iut) = 0<-^ulY,ul = M~iul) inn 

and 

u*{x) = ul{x) for all x € dQ. 
So by the comparison principle, for all i and for all e we have the upper bound 

Ui{x) > uj(x) for all x € 
Taking limits in e we can deduce that for all i 

(f^ii^) ^ Ui{x) for all X G dn. 

(b) Now consider as a lower barrier the solution of the problem: 

M~{uji) = inQ and coi{x) = (l)i{x) — ^'^4'j{x) indft 

So we have that for alii = 1, . . . ,d, 

•^"K - ^ «P < = M-{uji) in n 

and 

u^{x) — '^^Uj{x) = (jJi{x) for all x E dU. 

So by the comparison principle, for all i and for all e we have the lower bound 
ul{x) — 'Uj{x) > uJi{x) for all x G 

Taking the limit in e we can deduce that for all i 

Ui{x) — ^^■Uj(x) > ^i(x) — for all x € dil. 

Since by hypothesis we know that (f)i have disjoint supports, when 4>i{x) 7^ 0, 

Ui{x) > (t>i{x), 

and this proves the statement. 
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□ 



6. LiPSCHITZ REGULARITY FOR THE FREE BOUNDARY PROBLEM: PROOF OF THEOREM 12.51 

The regularity theorem of this section is a very important result relating the growth of one of the 
functions in terms of the distance of the function to the free boundary. In the proof we will need 
to use barriers, properties of subharmonic functions and the monotonicity formula introduced in 
[1]. In order to simplify the reading of the paper, the linear decay to the boundary is done in 
the beginning of this section and uses the construction of fundamental barriers that is done in 
Appendix E. The monotonicity formula is stated in Appendix [D] and the study of the decay 
for subharmonic functions supported in a small domains is presented in Appendix [Cj 

Lemma 6.1 (Linear decay to the boundary normalized). Let v be a non-negative continuous 
function defined in ft = Ba-{zo)\Bi{0), where cr < ^ and zq is, without loss of generality, a point 
on dBi{0) with = Cn, Cn the unit vector. Assume that, 

(1) M+{v) >{) infl, 

(2) v{x) < Ua in Q, 

(3) v{x) = on dBi{0). 

then, there exists a universal constant C, C = ^- — ,".a+i i such that, 

5~V 5 j 

v{x) < CUdist{x,dBi), 
when X & Sa with Sa ■= (0)\i?i(0)^ Pi {x = : \x' — Zq\ < |}. In particular, 

v{x) < CUdist{x,dBi), 

when X G Bz (zq). 



Proof. Consider the function v extended by zero to all i?o-(-Zo)- Obser ve that the extension still 
satisfies the hypotheses. We will take a barrier (p as in Lemma 



E.2 



with r = ^, ^ = I and 
M = C/|, that will be used as a model and will be sliding tangentially along dBi{0) in order to 
construct a wall of barriers (see Figure jTJ). With that purpose, take a family of balls {B^{y)}y 
such that y G B3a^{zQ) n dBi_sL{0). So the balls BsL{y) are tangent to dBi{0) and are inside -Bi(O), 

8 8 8 

where v is zero. For each ball consider (f> is such that: 

(a) (j){x) = Ua for x G dB5a_{y); 

8 

(b) 4>{x) = for x G dB^{y); 

(c) A^ + (0) <Oin B^^{y)\B^^{y); 

(d) g(x) = Ul^-^ when \x-y\ = f . 

5 \5> 



8 



Note that UyB5a_{y) C Bfj{zQ) and so for all y defined previously and for x G 9 \^UyB5a{y)j we 
know by hypothesis that 

v{x) <Ua. 
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We now apply the comparison principle for each barrier depending on y and respective ring 
B5a{y)\Bz{y), since v is a subsolution and <J) a supersolution for Ai^ , and we obtain that 

8 8 

<P{x)>v{x), x£dBs.{y)UdB^{y)^(l){x)>v{x), x £ Bs.{y)\B^{y). 
Hence, repeating this for all y we obtain that 

v{x) < (t){x), 



for all X €^ Scr = ^i^i+i (0)\-Bi(0) j D {x = {x',Xn) : \x' — Zq\ < |}. Taking in account that is 
radially concave, we also obtain that. 



(12) v{x) < U- ^ dist{x,dB,{0)). 



a 

For the final remark, observe that B^{zo) C Sa- □ 



Corolary 6.2 (Linear decay to the boundary). Let v be a non-negative continuous function defined 
in fifQ = Ba{tozo)\BtQ{0), where a < ^ <ii^d toZQ is, without loss of generality, a point on 5i?f(,(0) 
with = Cn, Cn the unit vector. Assume that, 

(1) M+iv) >0 in rito, 

(2) v{x) < Ua in Qt^, 

(3) v{x) = on dBtg{0). 

then, there exists a universal constant C, C = ^- — , such that, 

5~V 5 j 

v{x) < CiJdist{x,dBt,), 
when X £ with ■= (^B^^^a{0)\Btf^{0)^ fl {x = {x',Xn) ■ \x' — t^z'^l < |}. In particular, 

v{x) < CUdist(x,dBt,), 

when X £ Bs_{tQZQ). 



Proof. Consider the function 

u(x) = — v(xtn) 

Uto ' 

defined for x € = Ba(zo)\Bi{0), with a = f^. Observe that, u satisfies the hypotheses of 
Lemma |6. II with U = 1 and notice that cr < ^. Then, 

u{x) < Cdist{x,dBi), 

when X £ Sa- Substituting, we obtain 

^v{x to) < ^dist{tox, dBt,) ^ v{y) < CUdistiy,dBt,), 
Uto to 

for y £ Sa- □ 



We finally present the proof of Theorem 



2.5 
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Proof. (1) The proof is by contradiction. Let xo,zq S 0(supp ui) be two interior points to be 
characterized later. We will find a smooth function rj such that in Bs{zq) C i?ij(xo) touches the 
supersolution 

from below at zq, and simultaneously satisfies 

A^-(r?(zo)) >0, 

and this contradicts the definition of supersolution. In fact, we will find a positive universal constant 
Co such that if ui{y) = CR for some y G i?/j(xo) and C > Co then we have a contradiction. 

To simplify the notation let u = ui and v = '^^^lUk- Let us assume that the free boundary 
intersects the ball centered at the origin, 9(supp u) R -Bi (0) ^ 0. 

By contradiction, assume without loss of generality that u grows above any linear function in a ball 
centered at xq] more precisely, assume that for any constant M' and 

xq G 9(supp u) n Bi (0), 

2 

there exists y such that 

y £ Bii{xQ) and u{y) = M'R, 

where R < ^. Also, we may assume that R < 2d where d = dist(?/, Ssuppu) > (if not we can 
always pick another xq) and so we have 

y e Br{xo) and u{y) = M'2d = Md. 

For a later purpose we fix zq G dB^iy) n 9(supp u) (the closest point to y in the free boundary). 

As u > and u £ S*{X,A,0) we can apply Harnack for viscosity solutions on B^iy), 

sup u < c inf u 

and conclude that for any x G Bd (y) 

2 

(13) -Md<u(x)<cMd 

c 

Observe now that if uj(z) defined on Bd{y) is a solution to M~{uj) = then due to the invariance 
under translation by y (14), rotation by i? (15) and dilation by ^ and rescaling by d (16): 

(14) uj{x)=uj{x + y), X e Bd{0) ^ M-{oJ{x)) = M-{u}{z)), x e Bd{0), 



(15) uJ{x) = uj{Rx), x e Bd{0) ^ M-iuJix)) = M-{uj{z)), x £ Bd{0), 

(16) lJ{x) = ^oj{dx), X e Bi{0) ^ M-{uj{x)) = dM-{oj{z)), x € Bi{0), 

u){x) defined in -Bi(O) is still a solution to (uj) = with the direction e„, as we want. 

So we will prove this theorem using translation, rotation, dilation, and rescaling arguments on 
(ui, M2; ■ ■ ■ ) Ud)- In order to simplify the notation the new functions will always be denoted by the 
same name. 

Consider the functions u, v and u — v satisfying 

M-{u)>0, M-iv)>0, M-{iu-v)ix)) <0, 
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(see 2 in Lemma 5.1), and defined on an appropriate domain by translation (14), rotation (15) 



Cn and d is now 1. We wi^ 



dilation and rescaling (|16|) such that y is the new origin, the direction |^°_^| is now the direction 



1 still call zq the point in dBi{0) D dsuppu. By (13) we now have 

^^^Md<u(x) <\cMd, xeBi(0), 
dC d 2 



<u{x) <CM, xeBi{0). 



that is, 
(17) 

The rest of the proof consists of the following steps: 
First step: we will prove that for a certain M, a positive large constant, 

u{x)>Md{x,dBi{0)) X £ Bi{0)\Bi{0). 

^ V ' 2 

l-\x\ 

Second step: we will prove that, for a small p, 



(18) 



v{x)<^d{x,dBi{0)) xGSfi, 

\z\-l 



where Se := (0)\Si(0)) n {ix',Xn) G M" : \x' - z'q\ < f} (see Figure|7|. 




Figure 7. Barriers to control v. In the picture S is Sp defined in Lemma 



6.1 



Third step: Finally we will construct the function rj (see Figure [s]) to obtain the contradiction, since 



we will have for 5 < §: 



-M^{u-v)<0 for allx G ^^(zo); 

- r]{x) < {u — v) (x), for allx G Bs{zo) 

- v{zo) = {u-v) (zo); 
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- ?7 is a smooth function such that Ai^ {rj{zQ)) > 0. 
and by definition of supersolution this is impossible and the result follows. 
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Figure 8. Barrier function that touches u — v from below at zq. 



First step: By Lemma E.l (r = l,a = 1,6 = 2) there exist a subsolution ip of A4 on the ring 



Bi(0)\Bi(0) such that: 

2 

(a) ipix) = for G dBi{0); 

(b) i^(x) = for X G dBi(0); 

^ 2 

(c) M~{tp) > for x G Bi{0)\Bi{0); 

4 

(^) ^alr^ = ~ 2"°-i 17 ^ ^ c^-Si(O), where i/ is the outer normal direction. 



Due to the comparison principle applied in the ring Bi(0)\Bi(0) and by (17) one can conclude 

2 I 

that 

u(x)>i'(x), x£dBi(0)\JdBi(0)^u(x)>ij(x), x £ Bi(0)\Bi(0), 

2 2 

and also, since ^p is convex in the radial direction, and by (d), 

a M 

(19) > dist(x,5i?i(0)) X e Bi{0)\Bi{0). 

M 

Second Step: We do not know anything about the free boundary or the shape of the support of the 
other densities v = X^j-^i Ui. But we know that if zq is the point on the free boundary closest to 
then by the monotonicity formula 

dx ^ / — ^ — ^dx < 



since the norm of the densities is bounded, say by the fixed constant A^: 

C(n)||(n,z;)||t2(B^(0))<iV. 

Our goal in this step is to prove that v has to grow very slowly below a linear function with small 
slope, away from the free boundary. 

Let us consider the two possible cases around the point zo- 
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Second Step a): For a sequence of small radii, the measure of the intersection of the support 
of V with each ball is almost zero: there is a sequence of radii pk, Pk — ^ 0, such that 

\{vy^O}nBp,{zo)\<e\Bp,{zo)\, 

or, 

Second Step b): the measure of the support of v contained in a ball with any radius is not 
small, say, if for any radius p > 

\{v^O}nBp{zo)\>€\Bp{zo)\. 



The proof proceeds separately in each case but in both cases we will prove (18). 

Proof of Second Step a): If there exists a sequence of radius {pk)keN, Pk 0, such that 

\{v^O}nBp,{zo)\<e\Bp,{zo)\, 
we can consider a subsequence of {pk)k such that 

^ Pk ^ ^ Pk-l . / / 1 

■■■Pk+i < y < /Ofc < < Pfc-i • • • < Pi < 1- 

Since v is bounded in all domain, we have that there exists A^i such that 

sup v{x) < Nipi. 



Then, since v is subharmonic, by Proposition C.2[ we have that 

sup < eiVipi2". 



xeBpi (zo) 



Then, since 



M+iv) > M-iv) > ^M-iui) = 0, 



by Corollary 6.2 with If = eNi2'^^^, a = ^ and to = 1) we obtain that there exists a 
universal constant C : 

v{x) < C'eiVi2"+Mist(x,5Si), 
when X £ Spi with Sn := (0)\Si(0)) D {x = : \x' - Zq\ < Let e := eo 

2 2 \ 8 / 

be such that 



(20) 



eoC2"+^ < 



Then, for x G Sn 

2 



v{x) < —dist{x,dBi). 
Take as the next radii in the subsequence p2 such that Bp^{zQ) C Sn. Then, 



sup V{X) < —P2- 
xeBp,-,{zo) ^ 

So we can repeat the process, and so, again by Proposition 

iVi „ 

sup v[x) < eo— P2 2 . 

X&Bp2 (zo) ^ 



C.2 



with N2 = ) we have that 
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Then, by Corollary |6.2| with [/ = eo A^i 2", = f and to = 1, and by ([20j), 

v{x) < CeoA^i2"dist(x,aBi) < -^dist{x, dBi), 
when xG Sp2 with Sn := (B-^^p2{0)\Bi{0)) n{x = (x',x„) : \x' - z'q\ < f }. 

2 2 \ 8 / 

Again take as the next radii in the subsequence ps such that Bp^{zo) C Sp2. Then, 

sup V[X) < — P3. 
xeBp,^(zo) ^ 

Repeating a finite number of times this process, we have that, for x E Sn 

2 

Ni 

v{x) < ^dist(x,9Bi). 

Therefore, it is possible to find pi such that ^ < C=, and so 

v(x) <C=d(x,dBi(0)) xgSpi. 
M ^ V ' 2 

|x|-l 

In particular, 

v(x) <C=d(x,dBi(0)) xGBpi(zo). 
V ' « 

\x\-l 

Proof of Second Step b) : If for any p > as small as we want 

\{v^O}nBp{zo)\>€\Bp{zo)\. 

Since 

\{v / 0} n Bp{zo)\ = \{u = 0} n Bp{zo)\ > e \Bp{zo)\ , 
we can apply Poincare-Sobolev inequality to the function u. Also, by construction 

\{u ^ 0}\ = \{v = 0}\ > \B,{zo) n 5i(0)| > e |5p(zo)| , 

so we can also apply the Poincare-Sobolev inequality to v. Hence, 

v^{x)dx < C{n, e)p2 / \Vv{x)\^dx, 

Bp(zo) JBpizo) 



Then, 



u^(x)dx < C(n, €)p^ / |Vu(x)|^dx. 



[ t-2(x)dx<C(n,e)p2 /■ \Vv{x)\^ ^" _^ dx, 

JBpizo) JBpizo) \x-Zo\^ 

f u\x)d^<C{n,e)p' [ \yu{x)\ \ ^" ' , dx, 

JBJz^) JBJzn) k-2;or 



and 



v^ix)dx / u2(x)dx < C(n, e) p2"-4p8 ^_ 

B„(;.o) / XJBpizo) / 
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But we also know that u is controlled from below by the barrier function from the First Step, so if 
^ = 5i(0)\5i(0) 

/ V^(x)dx< / u2(x)dx< / u'^{x)(ix. 

JBp{zo)nA JBp{zo)nA JBp{zo) 

To estimate the integral of the barrier function we can use a linear function defined all over dBi (0) 
by the slope evaluated at zq. 

As in the first step, take M = 2<^-i % ^o be the absolute value of the slope of the barrier function. 
Note that M is arbitrarily big and so = will be arbitrarily small. Assume that the volume of the 

cone defined by the barrier function is bounded from below by cip'^p"'0Jn^ where ci depends on the 
area of the base of the cone, then, 

CiM p^p^-UJn . f .2 



< / V^(x)dx 

JB„(za)nA 



n JBp{zo)nA 

therefore, 



f 2/ ^ , W ClM^p2 "Wn\ . ^ 2n-4 Sat 

/ V {x)dx <C{n,e)p p N, 

JBpizo) J \ n 



so 



/ V {x)dx < ^2 = ^ ci, e)p ^ — ^. 

\JBp(zo) J CiM p'^p^OJn M 

Let y G Btp^^zq) such that B_p_{y) C Bp{zo). Due to the subharmonicity and positivity of v 

8 16 

v{y) < T f (x)dx < ( ) -r v{x)dx. 

JBpiy) \TqJ JBpiza) 

It) 



On the other hand, the Holder inequality yields 

2 



w(x)dx 

Bp{zo) 

All together, for any y G Btp{zq), we have. 



\ <i v{xfAyi. 
J JBp{zo) 



\y) < 16^" / v{xfd^ = le^"^— f v{xfdx, 
JBo(zo) P^^n JBo(zo) 



3p{zo) P"'^n JBp{zo) 

and hence 

1 fi^" AT AT 

v\y) < —C{a,n,c,,e)p^+^^ = le^" C(a:, n, ci, e) 

M ' V ' M 

4 

We conclude that v{y) < C2 = ^, for y G B7p{zo). Note that doing the same type of argument for 
r = pt, < t < 1 we can conclude that 

(22) y(y)<c2VN^\y-zo\. 
Since, v a subsolution of the positive extremal Pucci operator 

M+{v) > M-iv) > ^M-iui) = 0, 



A FREE BOUNDARY PROBLEM ARISING FROM SEGREGATION OF POPULATIONS 



37 



and we have by (22) that 



by Corollary 



6.2 



viy)<C2VN^^, x£Be{zo), 
with U = C2-^, a = ^, and to = 1, we obtain that 



v{y)<C^^distiy,dBi{0)), when y e S^. 



As before, Se is the portion of an annulus around -Bi(O) contained in Be{zq), 

Sa = (i?i+|(0)\i?i(0)) n {{x',Xn) G : \x' - z',\ < 
More explicitly, we have that for y G Bp(zq), 

(23) 



v{y) < C2VnI^- ^dist{y,dBi{0)) . 



5 1 _ fi^"+^ M- 

5 \5J 



12/1-1 



C 



Third step: Consider p the radii in Second Step a). Consider e to be eo in Second Step a), and C 

-8 a 



in (21) to be 



C = C2VN 



Putting together (19), (21) and (23), if A = Bi{0) D Bs{zo) and B = S D Bs{zo) 



u{x) > M(l - |x|) X £ A and 



C 

v(x) > — = (|x| — 1) X £ B. 



So letting rj be the radial symmetric smooth function defined by (see Figure [s]) 

r^(x) = (1 - |x|) + M(1 - \x\f, 

we can conclude that 

(u — v){x) > r]{x) X G Bs{zo) 

and that 

(n - v)izo) = rjizo). 
But notice that the Hessian of rj is given by 



Hmx) 



2M 







0^ 

-l+2M(r-l) 







-l+2M(r-l) 
r 








-l+2M(r-l) 



(where r = \x\) and so at zq because r = 1 we obtain M.~ {r]{zo)) = X2M — A{n — 1) > 0, since 
M is as big as we wish. But this contradicts the fact that n — w is a supersolution of A^~and this 
finishes the proof of (1). 

(2) Arguing in a similar way as in Proposition 4.13, in |2j, we should consider: 
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a) When2|x — y| <h, /i = max(dist(x, 9suppui), dist(y, Ssuppui)), y G i?h (x) so by interior 



regularity (Proposition 2.13 with / = 0) for properly scaled balls gives that 



\ui{x) - ui{y)\ < C ll-ui 11^00(5^(0)) 1^ - yh 

b) When 2\x — y\ > h note that for x in the support of ui and x is the closest point to x on 
the free boundary we have from the previous result that 

\ui{x) — ui{x)\ < < Ch. 

And so, by adding and subtracting ui{x) and ni(|7), where y is the closest point to y on 
the free boundary, and by triangular inequality, 

\ui{x) - ui{y)\ ^ \ui{x) - ui{x)\ ^ |m(x)-m(|7)| ^ \ui{y) - ui{y)\ 



\x — y\ \x — y\ \x — y\ \x — y\ 

^ 2|ni(x)| ^ 2\uiiy)\ ^ 
~ h h ~ 

Thus, 

\Wi\\Lip{Bi{xo)) = ll'"'illc">i(fii(a;o)) - ^' 
and the result follows. □ 
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Appendices 

Appendix A. Pucci Operators. General properties 

For the sake of completeness, we recall some properties of the Pucci operators in the following: 
M^(u)= inf aijDijiu) = inf TriAD'^u) = kS^ ei + a 

ei<0 ei>0 

where: 

^a,a is the set of symmetric matrices (n x n) with with eigenvalues in [A, A] for < A < A; 
ej eigenvalue of the matrix D^u\ 
Remark A.l. 

Observe that inf ag^;, ^ aijDij{u) < aijDij{u) < sup^g^^^^ aijDij{u). 
Let u, V be smooth functions and < A < 1 < A : 

- M-{u) < AAu < M+{u) ; 

- M-{-u) = -M+iu); 

- M^{u) + M^{v) < M^{u + v) < M^{u) + M^{v) and so M'^ is concave; 

- M^{u) + M^{v) < M^{u + v) < M^{u) + M^{v) and so is convex ; 

- < M^{u\) < AAu- < M.^{ul) =^ uj subharmonic in the viscosity sense ; 
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- < Xlt-^ (^i) — ■^{J2i^i) — ^ (X^i^i) ^ Z^j subharmoiiic in the viscosity sense; 



Appendix B. Fabes and Strook Inequality 
Here we state without proof part of Theorem 2 in [20]. For more details see also Appendix B in 

Lemma B.l. Let G{x,y) denote the Green's function for a linear operator L with measurable 
coefficients Lu = aij{x)DijU, [aij{x)] £ ^a,a- Then, there exist universal constants C and P such 
that whenever E C B^, and Br C Bi the following holds: 

2 

f G{x, y)dy <C I G(x, y)<iy for all x £ B^. 




Appendix C. Properties of subharmonic functions in thin domains. 

Lemma C.l decay for subharmonic functions supported in small domains). Let u be a non- 
negative subharmonic function in a domain that contains Bi{0). If for some small eo > 0, 

, \{u ^ 0} n Bi(0)\ 

sup -u < 1 and ,„ < eo 

Bi(o) l^iWI 

then, 

sup u < eo 2". 

Bi(0) 

3 



Proof. Let y be an arbitrary point in the ball Bi(0). Due to the subharmonicity and the fact that 

2 

u is a non-negative function 

1 



u{y) < 7 u{x)dx < 



Bi(y) 

2 



Bi{y) 



u(x)dx < 



Bi(0) 



i?i(o) n supp u\ 




Bi{y) 





sup u{x) 

Bi(0) 



and so by hypotheses, 



which gives the result. 



u[y) < jjyri — = eo 2 

(2) 



□ 



Proposition C.2. Let u be a non-negative subharmonic function in a domain that contains Bi{0). 
If for some p < 1 and for some constants N,eo >0 

\{u^o}nBp{ 



sup u{x) < N p and 



\Bpm 



< eo, 



then. 



sup n(j;) < iV/9eo2". 



x&Bp{0) 

2 
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Proof. Consider the function v defined on i?i(0) by 

The new function v satisfies 

sup v{x) < 1, 
xeBi{o) 



and so by Lemma C.l , we have that 



sup v{x) < eo 2". 
xeBi (0) 

Substituting v by its definition in terms of u gives that, 

1 

sup —u{px)<eo2"'4^ sup u{y)<Npeo2"- 
which gives the final result. □ 



Appendix D. Monotonicity formula 

Lemma D.l. Let u, v be two Holder continuous functions defined on Bi{0), nonnegative, subhar- 
monic when positive and with disjoints supports. Assume that belongs to the intersection of the 
boundary of their support. Let zq be the point in the free boundary close to 0. Then, the following 
quantity is increasing with the radius p and uniformly bounded, 



Jb,{zo) k - ^ol" ^ J \P'^ Jb,{zo) k - ^ol" 2 y 
since the norm of the densities is bounded, say by the fixed constant N ; 

Cin)\\iu,v)\\l,^s^^,)^<N. 

For the proof see page 214 on [5j and the proof of Lemma 7 (a) in [3]. 

Appendix E. Existence of Barriers 

Lemma E.l. Given constants < A < A and M,r, a,b, p > 0, ^ < r < p there exist a smooth 
function defined on Bp{0)\B^{0) and a constant c = c{a,a,b) and a universal constant a such 
that: 

(1) ^(x) = rMforxe dB^ (0) 

(2) ipix) = Oforx£ dBriO) 

(3) M-(ip) >Oforx£ Br(0)\Bij:(0) 

(4) ^ =cM when x G dBr{0), and c = -Ojj^^. 

Proof. Consider first r = 1 and M = 1, we will rescale afterwards. Let a, M2 > 0, a > n — 2 and 

ip{x) = Mi + M2- ^ 



x\ 



where Mi, M2 and a are such that the following conditions are satisfied: 
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(1) f{x) = when \x\ = 1 
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(2) ip{x) = 1 when |x| = ^ < 1 

(3) M-{(p) > in Bi{0)\Ba_(0) 
In detail, 

(1) (p{x) = when |x| = 1 ^ Mi = -M2 

(2) cp{x) = 1 when \x\ = ^ ^ M2 = = ^ > 
and so have that, 

(p{x) = -M2 + M2r^ with M2= " ^ . 

Note that if ^ is very small then M2 is very small too. On the other hand, the second derivatives 
of if are given by: 

dij(p{x) = — QM2 ^ij ^ a{—a — 2)M2XiXj . 

Evaluating the Hessian of at a point (r, 0, • • • , 0) one obtains: 

dijip = ij^j 

duif = M2a{a + l)r-"-2 

dii(p = -aM2r-"-2 i > 1 

And so by radial symmetry 

dijip{x) = i / j 
diMx) = M2 {a{a + 1)) |xr"-2 
dinp{x) = -aM2 i > 1 

which implies that the Pucci extremal operator is given by 

M-{(p{x)) = M2\a{a + 1) Ixp^-^ - A (n - 1) aMa Ixp""^ 
= M2a Ixr^^"^ (A(a + 1) - A (n - 1)) 
In order to satisfy (3) one needs that: 

A(n - 1) - A 
A 

which gives that 

,,1 . , a'" , A(n- 1) - A 

ipix) = -M2 + M2-r-[a with M2 = and a > — . 

\x\ 0" — a" A 

Notice that the normal derivative is: 

, , , ^ 1 

^(X) = -«M2^. 

And so when r = Ixl = 1, 



dip 
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Let C = —a-nr- 



Now let u s consider a dilation and obtain the result for general r and M = 1. Let <f{x) = r(p{^] 
Then (p defined on Br{0) satisfies: 

(1) <f{x) = when \x\ = r 

(2) <p{x) = r when \x\ 



a r 



(3) M-i<^)>0 in BriO)\Br^{0) 

(4) '^{x) = c when |a;| = r 

Finally for an arbitrary M, let il>{x) = Mip{x). Now is easy to check that the barrier function ^ 

satisfies: 

(1) ip{x) = when \x\ = r 



(2) ipi^x) = rM when |a;| 



ar 
b 



(3) > in Br{Q)\B9^{id) 

(4) |^(a;) = cM when \x\ =r 



□ 



Lemma E.2. Given constants < A < A and M,r,a,b, p > 0, < r < p there exist a smooth 
function defined on Bp{0)\B^{0) and a constant c = c{a,a,b) and a universal constant a such 
that: 

(1) ipix) = rM forxe dBr{<d) 

(2) Tpix) = Oforx€ dBriiO) 

(3) M+{ip) <Oforxe BriO)\B^^^{0) 

(4) ^ = cM when x G 55^o(0), where c = a- — y\a+i 



Proof. Consider first r = 1 and M = 1, we will rescale afterwards. Let a, M2 > a > n — 2 and 

ip{x) = Ml- M2j^. 

\x\ 

where Mi, M2 and a are such that the following conditions are satisfied: 

(1) (f{x) = when \x\ = | 

(2) ip{x) = 1 when = 1 



(3) M+if) < in Bi{0)\B^{0) 



A FREE BOUNDARY PROBLEM ARISING FROM SEGREGATION OF POPULATIONS 43 

In detail, (1) and (2) imply that 

a" 

Ml = 1 + M2 and M2 = 



b" - a" 
and so have that, 

1 a" 

ip{x) = I + M2 - M2r^ with M2 = - . 

|x| b°' — a" 

On the other hand, the second derivatives of (/? are given by: 

dijip{x) = aM2 6ij + a{—a — 2)M2XiXj \x\~°'~'^ . 

Evaluating the Hessian of at a point (r, 0, • • • , 0) one obtains: 

dij(p = i / j 

diiip = -M2a{a + l)r""~^ 

diiip = aM2r-"-2 i > 1 

And so by radial symmetry 

dij(p{x) = j 
dii(p{x) = —M2a{a + 1) 
dinp{x) = aM2 \x\~"~'^ i > 1 
which implies that the Pucci extremal operator is given by 

M+iifix)) = A (n - 1) aM2 \x\~'^~^ - XM2a{a + 1) jxl"""^ 
= Msa Ixp"-^ (A (n - 1) - A(a + 1)) 
In order to satisfy (3) one needs that: 

Mn — 1) — A 

A (n - 1) - A(a + 1) < a > 

A 

which gives that 

.X , , . 1 . , , . a° , A(n - 1) - A 

Gcjfa;) = I + M2-M2--T7 with M2 = ^ and a > — ^ ^ . 

\x\ b" — a°' A 



Notice that the normal derivative is: 



^-{x) = aM2 



And so when |x| = |, 

-^(x) = aM2 XT = XT- 

Let C = „_/°xa+l- 
6 \h) 

Now let us consider a dilation and obtain the result for general r and M = 1. Let <^(a;) = r</j(^). 
Then defined on i?r(0) satisfies: 

(1) (^(x) = when \x\ = r|, 

(2) (p{x) = r when |x| = r, 

(3) M+((^)<0in5,(0)\i?,|_(0), 

(4) ^{x) = c when |x| = r|. 
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Finally for an arbitrary M, let tpi^) = M(p{x). Now, it is easy to check that the barrier function ip 
satisfies 

(1) il){x) = when |x| = r|, 

(2) il^{x) = rM when \x\ = r, 

(3) M+iip) < in BriO)\BriiO), 

(4) = cM when |x| = rf . 

□ 
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